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K^ . Of late, the field of BFKL physics has been the subject of significant 

*vj ' developments. The calculation of the NLL terms was recently completed, 

Q>^ , and they turned out to be very large. Techniques have been proposed to re- 

■<!;;j- ■ sum these corrections. These lectures provide an introduction to the BFKL 

^^ I equation and some of the recent developments, using DGLAP evolution as 

the starting point. 

<3 ■ PACS numbers: 12.38.Cy 

Q^' 1. Introduction 

(— I I Some twenty-five years ago Balitsky, Fadin, Kuraev and Lipatov (BFKL) 

set out to determine the high-energy behaviour of the scattering of hadronic 
objects within perturbative QCD. They found terms going as (aglns)"', 
/\ ' where s is the squared centre-of-mass energy. Since Ins is large it can 

^ . compensate the smallness of Os and thus it was necessary to sum this whole 

series of leading logarithmic (LL) terms. The result was that the cross 
section should grow as a power of the squared centre-of-mass energy s [1]. 
For the values of as — 0.2 that are typically relevant, this power comes out 
as being of the order of 0.5. 

Over the past few years much experimental effort has been devoted 
towards observing this phenomenon, and the conclusion has consistently 
been that while the cross sections do rise, that rise is much slower than s^'^ 
(see for example [2-6]). 

The solution to this problem was to have been in the next-to-leading 
corrections to the BFKL equation, terms as(Qslns)", which have been 
calculated over the past ten years [7]. The various contributions were put 
together last year [8,9], and to the consternation of the community turned 
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out to be larger than the leading contribution, giving cross sections that 
were not even positive-definite [10,11]. 

These lectures will illustrate the origin of some of the main features 
of both the leading and next-to-leading BFKL equations, using as a basis 
the constraints provided by the DGLAP equation, and follow on with a 
discussion, based on [12, 13], of how to solve the problems that arise at 
next-to-leading order. 

After a brief definition of the problem in the next subsection, section 2 
discusses the DGLAP equation as relevant for high-energy scattering, and 
shows how it can naturally be extended to the give the BFKL equation [1]. 
This is followed by an illustration of the lack of agreement of the latter with 
experimental data. Section 3 derives the main features of the next-to-leading 
corrections to BFKL and discusses some of the problems that ensue from 
their inclusion. Section 4 looks at how one can go beyond next-to-leading 
order and section 5 concludes. 

1.1. The problem 

Let us first define a little more carefully ^2 

the problem to be addressed. We want to ^ ^ ^ 

study collisions of two perturbative hadr- / ^ 

onic objects, figure 1, where the squared 
centre-of-mass energy s is much larger than 
the typical transverse scales Q^, Qq of the 
two objects, which in turn are much larger 
than the QCD scale, A^, in order for the 

problem to be perturbative. This is of phe- 

nomenological relevance for certain features '~Z ^ 

of small-x deep-inelastic scattering (DIS) ^0 

at HERA, high-energy 7*7* scattering at Fig. 1. High-energy collision of 

LEP and the NLC, and configurations at two hadronic objects. 

the Tevatron and LHC involving jets that 

are widely separated in rapidity. It is also of theoretical interest since the 

large parton densities that arise at high energies can lead to novel effects 

such as parton recombination and multiple perturbative scatterings. 

2. Leading-logarithmic order 

2.1. Deep inelastic scattering 

Rather than entering straight into the problem of general high-energy 
scattering, it is helpful to consider first high-energy scattering in which 
one of the two hadronic objects is much smaller than the other, i.e. deep 



electron 





(a) (b) 

Fig. 2. (a) Deep inelastic scattering, (b) Cut ladder diagram for the evolution of 
the parton distributions. 



inelastic scattering, figure 2a. We have the collision of a proton (of mass M^ , 
equivalent to Qq of figure 1) with a photon of virtuality Q^ ^ Mp, which 
we will view as our second hadronic object. The photon-proton squared 
centre-of-mass energy is s. High-energy scattering in this system, s ^ Q^, is 
generally referred to as small-x scattering because Bjorken-x is = Q^/s -C 1. 

As is well known, to correctly treat such collisions it is necessary to 
resum terms (aslnQ^)"; because the smallness of ctg is compensated by 
the large size of InQ^. This is DGLAP [14] or collinear resummation, or 
renormalisation group evolution. 

The cross section is proportional to the quark distribution at scale Q^, 
which is related to the quark distribution at another scale Qq by 



xq{x, Q'^) = xq{x, Qo) + "s In — j 

'^0 
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In an appropriate gauge this can interpreted as the first in a set of ladder 
diagrams (figure 2b), whose rungs are strongly ordered in InQ^. 

To understand the type of ladder that dominates at small x, we need to 
look at the splitting functions. A quark ladder (with gluon rungs) involves 



iteration of the Pqq splitting function: 



P.(.)^^ 






(2) 



while a gluon ladder (with gluon rungs) involves the Pgg splitting function 



Pggi^) 



Ca 



vr 
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+ 5{l - z)Po . (3) 



At small z, Pqq is constant while Pgg grows as 1/z. So at small x, gluon 
ladders with repeated iterations of Pgg{z <^ 1) dominate, i.e. we have strong 
ordering in z. 

With this is mind, let us examine the properties of the unintegrated 
gluon distribution: 



Tix, Q2 



dg{x,Q'^ 
' dQ2 



(4) 



and start with a simple (though not entirely physical) initial condition 

Q2^W(x,g2) = e(i_a,)e(^-i^ , (5) 

where the Q^ factor is included on dimensional grounds {g{x^ Q"^) is dimen- 
sionless). Using the purely gluonic DGLAP equation in differential form, 



Q^x 



dg(x,Q2 
dQ2 



as 
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z \z 
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and rewriting it in terms of the unintegrated gluon distribution, we obtain 
the first-order contribution to J^, 

where as = asC^/vr has been introduced as a notational shorthand and a 
factor Q{Q^ — Qg) is implicitly understood to be contained in the result. 
We retain only the 1/ z part of the splitting function because the other parts 
lead to contributions lacking the factor In \/x and so much smaller than (7). 



The second-order contribution is 

"^ dz2 



Q2^(^)(x,Q2 



Os 
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^m^lin^^ 
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By iteration one sees that the O (a") contribution is 



[niy y^" "' X "' Qf 



(8) 



(9) 



Since every power of as is accompanied by two logarithms, this is referred 
to as a double-logarithmic (DL) series. It resums ladders in which there is 
strong ordering of both the transverse and longitudinal momenta along the 
ladders: /c^/A;,^_]^ ^ 1 and Zn <^ 1 respectively. 

2.2. Summing the DL series 




Fig. 3. Graphical depiction of the integral equation (11) 



Our DL series happens to be related to the series for the modified Iq 
Bessel function [15]. Using the asymptotic expansion for /q gives us the 
result that 



Q^T{x,Q^)^exp 



iln^' 
X Q'q 



2\ /as In — In -^ 
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However it is useful to develop a more general method of summation, one 
which will be applicable also later on. Accordingly, we formulate the prob- 
lem as an integral equation 



Mz /-^'dP 



T{x,Q')=T^'\x,Q') + a,l^^l ^.^g. 



t 



(11) 



which is depicted graphicahy in figure 3. It can be diagonahsed by taking 
Melhn transforms with respect to both x and Q"^, 



2\ 7 



■^(a:,Q')-y 2vr^^" i 2vr^Q2(QgJ -^7,- > 
with the contours running parahel to the imaginary axis, giving 
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The pole in 7 is conjugate to the DGLAP logarithm of Q^ and the pole in 
(J conjugate to the logarithm of x. Eq. (13) is easily solved: 
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oj!F'- 
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With the initial condition (5), we have J--f,uj = l/i^T- The inverse Mellin 
transform with respect to lo is carried out by closing the uj contour to the 
left in eq. (12), leaving us with 



d7 _^ 



«^(-«' = ./2i^"^l% 
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The integrand has a saddle-point at 



as In 1/x 



and a saddle-point evaluation of the integral gives 

1/4 



(16) 



Q^J^ix,Q^)^ 
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This result was first obtained twenty-five years ago by De Rujula et al. [16]. 
Its main feature is that the gluon distribution rises at small x, with an 
effective power 



"^eflf 



' Og In 1/x 
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Fig. 4. The effective power of the rise of F2 [2]; x and Q^ values are correlated so 
that higher Q^ also means higher x. 



which decreases as one moves towards smaller x values, and increases to- 
wards higher Q^. The HI data for uJcS shown in figure 4 illustrate precisely 
this trend. Detailed comparisons of DGLAP-induced rises of F2, including 
full splitting functions and a running coupling, such as those performed by 
Gliick, Reya and Vogt [17] and Ball and Forte [18], show remarkably good 
agreement with nearly all the small-x F2 structure function data, even down 
to Q2 ^ 4 GeV2. 



2.3. BFKL 

The above arguments are relevant in a limit where both 1/x and Q'^ IQ\ 
are large, i.e. when we have strong ordering in both longitudinal and trans- 
verse momenta. But when the ends of the chain have similar transverse 
momenta, Q^ ~ Qq, there is no longer any reason for transverse momenta 
along the chain to be ordered. Double logs no longer dominate the cross 
section and we have to sum all leading (single) logarithms (LL) of x. 



Og In ■ 



1 
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This is done by the BFKL equation [1], which can be derived in a number of 
ways. Since these notes are intended only as an introduction to the BFKL 



equation, rather than engaging in its derivation we will try to deduce its 
main characteristics from simple physical arguments.^ 

In the previous section we had the following integral equation for the 
gluon density (11): 

T{x, Q2) = ^(0) (^, q2) + I^I dk'K{Q^ k^) T (^, fc2) , (20) 

where the DGLAP kernel K was just 

K{Q^, k^) = %<d{Q^ - k^) , valid for Q^ y^ k\ (21) 

Since the BFKL kernel should be valid for any ratio of transverse scales 
it must have the same limit for Q^ ^ fc^, and additionally correctly treat 
situations in which Q^ is of the same order as, or much smaller than k'^. We 
can deduce its form in the limit k'^ S> Q^ by the following argument. The 
scattering of a big object off a small one, or of a small object off a big one, 
must have the same cross section, and both situations must be correctly 
described by the BFKL resummation. Therefore the BFKL kernel must 
be symmetric under the interchange of Q^ and /c^. So when k"^ ^ Q^ (the 
anti-DGLAP, or anti-collinear limit) we have 

K(Q^, k^) = ^, valid for k'^ :^ Q^ . (22) 

k'' 

If we approximate the full kernel just by its collinear and anti-collinear 
limits, then we have 

Following the treatment of the previous section, we will need its Mellin 
transform, 

X^°"(7) = - + T^, (24) 

7 1-7 

where, by convention, the leading factor of as has been left out; xil) is 
usually referred to as the characteristic function of the system. The I/7 
term was present also in the pure DGLAP case, and comes from the collinear 
limit. The 1/(1 — 7) term comes from the anti-collinear limit. The symmetry 

^ For the interested reader one of the simplest fuU derivations is perhaps to be found 
within the dipole formulation [19]. A wide ranging introduction and discussion of 
many aspects of BFKL physics can be found in [20] . 



7 <-> 1 — 7 is a direct consequence of the symmetry under the exchange of 
the two transverse scales. This can be seen expHcitly from the definition of 
the Melhn transform, eq. (12), where 






1 



Ql 



(25) 



What we have neglected in our collinear + anti-collinear approximation 
is the correct treatment of the kernel for k of the same order as Q. This is 
given by the full BFKL kernel [1] (with integration measure cPk/ir because 
the azimuthal integration now matters): 



K(Q2,fc2) = a, 



1 



|g-fc|2 



<5(Q2 




Its Mellin transform (again leaving out the overall factor of Og) is 

x(7)=2^(l)-^(7)-V'(l-7), 



(26) 



(27) 



where 'ip{x) = dlnr(x)/dx. Noting that — '(/'(t) = I/t + C'(l) for small 
X, we see that the full xil) h^-s the same polar structure around 7 = 
and 7 = 1 as our approximation (24) reflecting the fact that the collinear 
and anti-collinear limits are the same. The two characteristic functions are 
shown in fig. 5, which illustrates their very similar shapes: they differ by 
little more than a constant. 




Fig. 5. The full and collinear ly-approximated characteristic functions. 
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Fig. 6. The real part of the LL BFKL characteristic function in the complex plane 
and the integration contour for the inverse Mellin transform. 



The procedure for obtaining the BFKL cross section is analogous to that 
used in the DGLAP case, with I/7 replaced by xil)- We start with the 
Mellin-transfornied integral equation 



^7,- = ^i'l + -X(7)-^' 



7,tj ) 



and solve for J^. 



J,U1- 



T. 



ujT^ 



7,0; 



7,0; 



(28) 



(29) 



uj - asx(7) 

The inverse Mellin transform with respect to uj is once again trivial and 
gives 



T{x, Q^ 



d7 

27ri 



exp 



asX(7)ln- +7 In —2 



ojT. 



(0) 



7,0; 



w = OisXil) ■ 
(30) 



The real part of the characteristic function x{l)^ together with the integra- 
tion contour are shown in the complex plane in figure 6. For large In 1/x, 
the integrand in (30) is dominated by behaviour of x ^-nd has a saddle-point 
close to 7 = 1/2, which causes the gluon distribution to grow as 

™-asx(|) 1 



T{x,Q' 



27rasX"(^)lni ^^^ ' 



(31) 
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where % is the second derivative of x with respect to 7. Thus the gluon 
distribution (and in general, high-energy cross sections) should grow as a 
power of X determined by the minimum value of x{l)^ which is x(l/2) = 
4 In 2. For the phenomenologically reasonable value of Qg ~ Os = 0.2 this 
gives a power of about 0.5. 

It suffices to look back at figure 4 to see that this is incompatible with the 
rise seen in the bulk of the structure function data. Some care is needed in 
interpreting this disagreement: in considering the structure function data, 
we are trying to apply perturbative QCD to a problem which is inherently 
non-perturbative (the scale Qq does not satisfy our requirement Qq ^ A^). 
However BFKL also predicts scaling violations of the F2 structure function 
[21], and this prediction can be shown not to depend on the properties 
of the non-perturbative region [13,22]. Essentially, regardless of the input 
distribution, the scaling violations quickly lead to a structure function which 
rises with a power 41n2as and so is incompatible with the data [23]. 

There exist other, theoretically cleaner, tests of BFKL. Generally they 
involve selecting a process with two hard hadronic probes, such as jets or 
a virtual photon, separated by a large rapidity (or equivalently having a 
large centre-of-mass energy). The requirement that both probes be hard 
ensures that one can reasonably apply perturbation theory^ (unfortunately 
it generally also makes the experimental measurement much harder). A nice 
example of such a test is the collision of two virtual photons as measured 
recently by the L3 [4] and OPAL [5] collaborations. The L3 data are shown 
in figure 7. The data are significantly higher than the one-gluon estimate 
(i.e. the prediction without BFKL resummation) . On the other hand the 
LL BFKL predictions clearly overshoot the data. The L3 collaboration 
perform a fit to the data in order to determine the power of the high-energy 
growth, and quote a preliminary result of 0.29 it 0.025 for scales in the range 
of 3.5 to 14.5 GeV^ [24]. ^ 

The same conclusion of incompatibility with LL BFKL comes out from 



Though not strictly the subject of this presentation, an exposition of BFKL physics 
would be incomplete without at least some mention of diffusion. Because transverse 
momenta are not ordered, small- a; evolution leads to a rand om w alk in Infct- The 
mean width of this random walk — diffusion — increases as Vms, and at very large 
s eventually enters into the non-perturbative region. Thus, no matter how large the 
transverse scales of the scattering objects, there is always an energy beyond which 
perturbation theory loses its predictive power. 

This result should probably be interpreted with some caution because the formula 
used to carry out the fit assumes the LL normalisation with four light-quark flavours 
(whereas both the NLL corrections and the charm mass probably have a significant 
effect on the normalisation). A fit leaving the normalisation as a free parameter 
leads to a similar power but with a much larger error. One should bear in mind that 
because of the 'limited' energies at LEP, the Q^ values (between 3.5 and 14.5 GeV'^) 
are probably on the border of the region that can be considered perturbative. 
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Fig. 7. The cross section for 7*7* collisions as measured by the L3 collaboration 
[4,24]. The mean Q^ values for the three energies are 3.5, 14 and 14.5 GeV^ 
respectively. 



looking at the interactions between a jet and a virtual photon in DIS [6], 
a measurement referred to as the forward-jet cross section, because of the 
position of the jet in the detector. 



3. Next-to-leading corrections 

All along, while the various experimental tests of LL BFKL were being 
carried out and refined, the calculation of the next-to-leading logarithmic 
corrections to BFKL was in progress. The next-to-leading terms are those 
suppressed by a power of as relative to the LL series: 



aslaslns)"" . 



(32) 
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In terms of the notation developed so far, this corresponds to working out 
the NLL corrections to the characteristic function X; ^-e. finding xi, where 

asX(7) = asXo(7) + ajxih) + O (a^) . (33) 

The determination of xi took close to ten years [7], and was completed quite 
recently [8,9]. 

Rather than trying to reproduce parts of that derivation, we will adopt 
the same approach that was used in the previous section, namely to de- 
duce the structure of the NLL corrections through a study of the collinear 
limit and symmetry requirements. This will translate to determining the 
divergences around 7 = and 7 = 1. 

We will examine three main contributions: those from the running cou- 
pling, the non-singular (at small z) part of the splitting functions and the 
choice of energy scale. 

3.1. Running coupling 
The QCD coupling runs as 

where b = 11/12 — nj/6. What sort of higher order contribution will this 
lead to? The DGLAP equations tell us that in the right-hand graph of 
fig. 3, when Q^ ^ k'^ the correct scale for the coupling is Q^. By symmetry, 
when k'^ ^ Q^, the correct scale is /c^ — i.e. in the collinear limit the correct 
scale is the larger of the two scales involved. So our collinear approximation 
for the kernel, eq. (23), becomes 

K [Q,k)=as[Q) ^2 \-o:s{k) p • (35) 

The Mellin transform of the first term just gives as{Q'^)/'y, as before. For 
the second term, we re-express as{k'^) in terms of as((5^) in order to extract 
a factor of as{Q'^) in front of the whole result. Expanding to second order, 
and taking the Mellin transform, gives 

which is just the anti-collinear part of our LL result plus a running-coupling 
NLL contribution 

^5 = -(r^- (^^' 
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The lack of symmetry 7 <-> 1— 7 is due to our choice to extract an asymmetric 
factor of as{Q^) in front of the answer. 

What is the uncertainty on our cohinear approximation for Xi? The 
scheme of a^ is not defined, corresponding to an uncertainty proportional 
to xo- Nor do we a priori know the correct scale for branchings when fc^ 
and Q^ are of the same order. So the overall uncertainty is a function with 
at most single poles at 7 = and 7 = 1. 




«s/t^7 



asAi 



Fig. 8. A sequence of a 
small- a; and a non-sinall- 
X branching. 



3.2. Splitting function 

In section 2 we used only the part of the gluon 
splitting function that is singular at small z. At 
NLL, we need to include the full splitting function 
(3). Its Mellin transform (with respect to x) is 



99 



-+Ai(w), 



where (for Uf = 0) 



Aiiiv) 
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(38) 



(39) 



To get the NLL correction we consider a sequence 
of two collinear branchings, fig. 8, where one of 
the branchings is a small-x branching, giving a factor Os/ujj and the other 
is a non-small-x branching, giving a factor 03^1/7. Remembering that 
convolutions in x and fc^ space translate to products in the uj, 7 Mellin 
transform space, our integral equation (28) receives a contribution 



as "8-4 



U! 7 



2 •' 7,t<J • 



(40) 



There is a corresponding term for a pair of anti-collinear branchings, so that 
the splitting-function contribution to xi is 



xr(7) 



7^(1 



Ai 



1?' 



(41) 



where Ai = ^i(O) = —11/12. Actually this is only the nj-independent part. 
For non-zero n/ there are contributions coming from the n/-dependence of 
Pgg and from diagrams involving the convolution of Pgq and Pqg . 

As in the running coupling case we have an uncertainty on this result, 
which can arise for example from the combination of a collinear and an 
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anti-collinear branching, and thus is once again at the level of a function 
with at most single poles at 7 = and 7 = 1. 

3.3. Energy scale terms 

A more subtle source of NLL corrections comes from the so-called 
energy-scale terms. At leading order one resums terms 

asln^") , (42) 

where so can be chosen arbitrarily. Changing so is equivalent to introducing 
a whole set of higher order terms. For a symmetric treatment (with respect 
to Q and Qq), a natural choice is sq = QoQ. Let us then consider what 
happens when Q ^ Qq. As we obtained in section 2.1, the leading terms 
are 

1 /_ s Q^N" 

"n\2 "s In — — In -2 . (43) 

But for DIS-like situations, Q ^ Qo) we usually express our results as a 
function of x = Q^/s and Q'^ IQ\. Let us do that for the n = 2 term: 

Uasln-^ln^) =\[a,\u-\^^ + 1 «2 ^^ i ^^3 g + NNLL . 
4 V <5<5o Ql) 4 V X Ql) 4 ' X Ql 

(44) 

Of particular interest is the second term on the RHS because it has more 
collinear logs than powers of a^ (it contains a double collinear logarithm for 
a single power of Og). But we know from renormalisation group constraints 
that the cross section written as a function of x and Q^ /Q^ contains at 
most as many collinear logs as powers of a^. Therefore for the result to be 
consistent with the renormalisation group, the next-to- leading corrections 
must be such as to cancel the second term on the RHS of (44), i.e. they 
must contain a term 

-^a^lnsln^g. (45) 

4 (^0 

In Mellin transform space this will correspond to a contribution to xi which 
is proportional to 1/7^. To obtain its coefficient it is not sufficient just to 
take the Mellin transform of (45), because not all of the correction expo- 
nentiates (i.e. should be incorporated into x) — foi' example some of it is 
to be associated with the initial condition. 
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Instead, now that we know what kind of answer to expect, let us directly 
consider the problem in Mellin-transform space. We start with a result 
written for energy scale QQo {cf. eqs. (12) and (29)), 





" 1 fQ^y uj^\% 


Q2 [q2) u;-a,xil) 


and then note that since 




[qQoJ Wo) 


W) UgJ ' 



(46) 



(47) 



rewriting (46) with energy scale sq = Q^ is just equivalent to shifting 7 
7- f inx and ^■(O): 

dc^ [ d-i ( sY I (Q^y ^-^7%,- 



If we now expand x(7 ~ ^) iii powers of ag, recursively using the relation 

uj = OsX we get 

X(7 - |) = X(7) -^ + {al) . (49) 

In the coUinear limit (7^0), since xil) goes as I/7, the O (ag) piece has 
the behaviour 

This is the analogue of the In Q^ term seen earlier (44) and it must be 
subtracted from x at scale sq = QQo in order for the collinear limit with 
energy scale Q^ to be free of unwanted double collinear logs. There is an 
analogous 1/2(1 — 7)'^ piece to be subtracted for the anti-collinear limit to 
be correct (i.e. free of double anti-collinear logs for energy scale sq = Qq). 
Overall therefore we have the following NLL energy-scale corrections (for 
So = QQo)- 

^1° = -^- 2(1-^)3- (^1) 

As was the case for the running coupling and splitting function terms, our 
analysis leaves us with an uncertainty which amounts to a function with at 
most single poles at 7 = and 7 = 1. 
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3.4- Putting things together 

Putting together eqs. (37), (41) and (51) gives us the fohowing answer 
for the colhnearly-enhanced part of the NLL corrections (nj = 0): 



coll 



Ai Ai - 6 1 1 



Xi 



^^> 72 ^ (1 - 7)2 273 2(1 - 7)3 ' ^ ' 



with Ai = —11/12. Our cohinear approximation guarantees the correctness 
of the coefficients of the cubic and quadratic divergences at 7 = and 7=1. 
The true NLL corrections as assembled in [8, 9] are, in the ms scheme 
and for ny = (the Uf dependence turns out to be small) 

7r2cos(7r7) / 2 + 87(1-7) 

Xi[l) = - , . 2/ — TV, ^TT -^ + 



4sin2(7r7)(l - 27) V (3 - 27)(1 + 27) 

-5{x§(7)-*'(7)+*'(l-7)) + ^ + ^^^ 

'67 7r2\ _ 3 ._ 7r3 



+ l36-l2J»<^' + 2';P' + 4ilK^-*<^'' <^=" 



where 

<^(7) = E(-i)" 



n=0 



V;(n + l + 7)-V^(l) V>(n + 2-7) - V'(l) 
(n + 7)2 (n + 1 — 7)2 



(54) 



It is possible to make a direct identification between parts of (53) and (52) 
in terms of the coefficients of the double and triple poles. The first line 
of (53) is identifiable with the Ai piece of (52) and so originates from the 
splitting function. The running-coupling dependence enters through the 
first term on the second line of (53), while the energy-scale dependent piece 
is formed by the last two terms of that line. The remaining terms are free 
of double and triple poles. Of these terms, so far only the first one on 
the third line of (53) has been understood: it is associated with the fact 
that the natural scheme for processes involving soft gluons is the CMW 
or gluon-bremsstrahlung scheme [25]. When writing an answer in the ms 
scheme this leads to a correction term which is (67/36 — 7r2/12)as times the 
leading order result. 

Figure 9 shows the full xi together with our collinear approximation. 
There is a remarkable similarity between them: in the range < 7 < 1 they 
never differ by more than 7%. Possible reasons for the surprisingly good 
agreement will be discussed later, in section 4. 
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Fig. 9. A comparison of our coUinear approximation and the full result for xi; 
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Fig. 10. The real part of the LL + NLL characteristic function for as 
together with the integration contour, in the complex plane. 



0.2, 



3.5. Consequences of the NLL corrections 

Figure 9 shows that the NLL corrections to x are very large. For 7 = 1/2 
we have (again for Uf = 0) 



x(l/2)=xo(l/2)(l-6.47as) 



(55) 
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Thus for as = 0.2, the predicted power is negative, at about —0.16, which 
is no more in agreement with the data than the leading power. 

What is even more worrisome is that the structure of the characteristic 
function changes radically. In figure 6, there is a single saddle-point at 
7 = 1/2, i.e. on the real axis. With the NLL contributions included, there 
are now two saddle-points, at complex values of 7 which we will call 7, 7* 
[10]. Since the cross section goes as 



«sX(7) I /g2\7 

QQoJ Q^ [qI 



^is,Q^Ql)-{7^] 7^(^) +7-7*, (56) 



the fact that 7 is complex means that the cross section oscillates as a 
function of IuQ'^/Qq (it remains real because of the contribution from the 
complex-conjugate saddle point). This behaviour occurs as long as x has a 
negative second derivative at 7 = 1/2, which is the case for a^ ^ 0.05. In 
other words there exists no phenomenologically accessible domain in which 
the inclusion of the NLL corrections gives a sensible result. 

4. Beyond NLL 

The solution is bound to lie with higher orders. Shortly after preliminary 
results on xi had appeared, it was suggested that stable predictions might 
be obtained by inclusion of the NNLL and NNNLL terms [26]. But remem- 
bering that the LL calculation took about a year, and the NLL calculation 
ten years, a reasonable estimate for the time to calculate the NNLL terms 
might lie somewhere between an arithmetic (19 years) and a geometric (100 
years) extrapolation. Even were these contributions to be calculated, there 
is actually no guarantee that the resulting series would converge for the 
values of as of interest! 

So the only option left is to try and guess the higher-order terms and 
then to resum them (we are now talking about the resummation of a re- 
summation) . The question is whether there is some reliable way of guessing 
them. Various approaches have been investigated [12,13,27-29]. I here will 
advocate a method closely related to that used in the previous section to 
estimate the NLL corrections — namely a method based on the study of 
the collinear limit [12,13]. 

We have already seen in the previous section that a study of the collinear 
limit is a powerful tool. For the NLL characteristic function it gave us the 
cubic and quadratic divergences at 7 = and 7 = 1, and in the range 
< 7 < 1 reproduced the full answer to remarkably good accuracy. How 
come? 

There is a temptation to argue that since 7 = 1/2 is a moderately small 
number, one can legitimately carry out an expansion in powers of 7 and 1 — 7 
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(including the first two terms in the expansion for xoi i-^- the two poles and 
a constant, also seems to do quite well, reproducing the full answer to within 
about 8%). A slightly better motivated argument might be the following. 
For 7 = 1/2, the transverse momentum integrals in the Mellin transform 
converge quite rapidly and so one might not expect a collinear approxima- 
tion to work too well. However at higher orders, pieces of the transverse 
momentum integrals are accompanied by logarithms of transverse momen- 
tum. These have the effect of shifting the dominant part of the integral 
out towards more collinear regions, where the collinear approximation itself 
becomes better. 

So there are reasons to believe that coUinearly-enhanced contributions 
might give a significant part of the higher-order corrections even beyond 
NLL.*^ A more general justification for carrying out the collinear resum- 
mation is that one wants to be able to use one's answer in the collinear 
limit. Since the collinear limit involves taking 7 close to zero, where higher 
orders involve successively more divergent terms, the only hope of a sensible 
answer there is a collinear resummation. 

The determination of the collinearly enhanced corrections can be divided 
into two parts. The first deals with terms in the same class as the 1/7^ terms 
in xi, single collinear logarithms originating from the splitting function and 
running coupling; the second addresses the terms in the same class as the 
1/7^ term in the NLL result, namely double collinear logarithms. We will 
consider only an outline of the method. The interested reader is referred to 
the original references [12, 13] for the full details. 

4.1. Single collinear logs — running coupling and splitting function terms 

It is fairly straightforward to calculate the collinear N"LL contributions 
to the BFKL kernel from splitting function and running coupling effects. 
One just takes diagrams such as figure 8 with an arbitrary number of non- 
small-x emissions, inserting and expanding the appropriate running coupling 
for each branching (the answer is given in [13]). One is then left with the 
tricky problem of resumming the resulting set of terms. 

An equivalent approach essentially treats the small-x and non-small-x 
branchings on a more similar footing [13]. In eq. (28) we have a factor as/to 
coming from the 1/z part of the Pgg splitting function, and a factor of xil) 
from the transverse structure of the branching. In the collinear limit, we can 
replace l/to with the full splitting function, P^ . So the collinear behaviour 

* This statement should really be restricted to those corrections that can be associ- 
ated with a single ladder (referred to as i-channel iteration). Actually at NNLL, 
corrections arise associated with the presence of two ladders (the start of s-channel 
iteration), i.e. saturation, or unitarity corrections. Our aim here is to understand the 
high-energy behaviour of a single ladder. 
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of the as,x/^ factor becomes (for DIS energy scale, sq = Q'^) 

^sP^g as{Q^)l+ujAi{Lo) 1 + ujAi{lo) ,^. 

— = => X^ , 57 

7 cj 7 7 

where we have used eq. (38) for P^ . The correct scale in the branching 
is Q^ (the largest scale in the problem) so there are no running coupling 
corrections. 

The situation in the anti-collinear limit is very similar except for an issue 
related to the running coupling: the appropriate scale for the branching is 
not Q^, but k"^ (referred to fig. 3). If we want to extract a factor of as{Q'^) 
the difference of scales must be taken into account. It turns out [13] that 
this can be done at all orders by replacing Ai with Ai — b, so that the 
resummed anti-collinear structure is (with the anti-DIS energy scale choice, 
so = Ql) 

«s(Q^) 1 + ^(^1 (^)-fc) ^ ^^ l + u{Ai{u;)-b) ^ ^^g^ 

u 1 — 7 1 — 7 

The reader can verify that substituting to = OsXo into the expressions for x 
in (57) and (58) reproduces the correct NLL collinear-enhanced terms. 

4-2. Double collinear logs — energy-scale terms 

We have just given resummed answers for the collinear and anti-collinear 
behaviours of the kernel with DIS and anti-DIS energy scale choices respec- 
tively. We really want the answer for a common energy scale, say sq = QQq. 
We saw in section 3.3, that changes in sq introduce higher-order double- 
coUinear logs (the 1/7^ and 1/(1 — 7)^ terms). The higher-order corrections 
had to be such that for energy scale sq = Q^ there were no such terms 
around 7 = 0, and similarly around 7 = 1 for energy scale sq = Qq. For 
•So = QQoj appropriate double-collinear log counterterms had to be included 
in order to satisfy the conditions for the other energy scales. One can work 
out, order by order, the counterterms for higher kernels, but it soon gets 
tedious. In any case one finds that the resulting series of terms is divergent 
for reasonable values of Qs- 

The solution [12] exploits the fact that a change of energy scale corre- 
sponds to a shift of 7 by an amount proportional to u (cf. section 3.3). For 
energy scale sq = QQo, one writes a leading-order kernel with the following 
structure 

X^ = 2V'(l)-V(7+|)-V'(l-7 + |), so = QQo, (59) 
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originally discussed in [30]. Changing energy scale to sq = Q^ corresponds 
to the shift 7^7 — ^ (c/. section 3.3), and we have 

X^ = 2^(l)-^(7)-^(l-7 + c^), so = Q\ (60) 

Remembering that — V'(7) — I/7 for small 7, and iteratively solving for 
cu = OisX as we did in section 3.3, we find an answer which is free of sin- 
gularities stronger than I/7, and so free of spurious double collinear logs. 
The procedure can be repeated for energy scale sq = Qq, expanding around 
7 = 1 and one finds an answer free of spurious double anti-collinear logs. 
Expanding (59) to order Os gives exactly the same triple poles as in (51). 

4.3. The full resummed answer 

Let us first see how to correctly include the energy-scale resummation 
in the full kernel. We start with the modified LL characteristic function, 
Xq, eq. (59) which, as we have just seen, is free of unwanted double (anti) 
collinear logs for the (anti) DIS energy scale choice; Xo contains NLL cor- 
rections, 

^(-V'(7)-^'(l-7)), (61) 

which must be subtracted from xi to avoid double counting: 

Xi = Xi-f (-V''(7)-^'(l-7)). (62) 

The quantity xi still has quadratic and single divergences at 7 = 0, 1. 
In analogy with the single divergences in xO: these need to be 'shifted' in 
order to avoid spurious double-collinear logs at higher orders when changing 
energy scale. This is accomplished by subtracting unshifted divergences and 
replacing them with shifted divergences: 

Xr = XI - ^i(0)V'(7) + ^i(^)V'' (7 + I) - (^i(O) - b)^P'{l - 7) 

+ (^i(^) - b)i;' (1 - 7 + I) + !^ (;t- - ^,) . (63) 

Here we have chosen to use ^'(7) and —ipij) (in xo) a-s our quadratic and 
single 'divergences to be shifted'. We could equally well have used 1/7^ 
and 1/7 respectively. The difference in the final result would amount to 
collinearly suppressed NNLL terms. The reason for including ^1(0;) in the 
shifted poles is discussed shortly. 
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To resuni the splitting- function and running coupling effects, we have to 
ensure that x has the following structure around 7 = and 7=1, 

xIt,'^)- r^, 7<1, (64a) 

7+2 

1 + uj{Ai -b) , , , 

X{l,u^)^ , 1^ , 1-7«1, (64b) 

1-7+2 

where the poles have been shifted compared to eqs. (57) and (58) to take 
into account that they have been written for energy scale sq = QQq. This 
can be obtained by writing 

x(7,u;)=Xo+^4- (65) 

Xo 

Since to = asXo + ^ ("^s) ^^^ expansion of x to order as is correct. Addi- 
tionally the ratio xf/Xo contains (shifted) single poles at 7 = and 7 = 1 
with coefficients Ai and Ai — b respectively, as required by eqs. (64). The 
full uo dependence of Ai is included through the Ai{u;) factors eq. (63). 

A point to note is that (65) is no longer an expansion in as, but rather in 
uj. For this reason this resummation technique is known as the w-expansion 
[13]. 

4-4- Results 

Figure 11 shows various BFKL exponents as a function of as, including 
the LL and NLL results for reference. The quantity labelled ujg is the 
minimum value of lj = asx(7)'^)) ^^^ ^^ such corresponds to the exponent 
expected for the gluon Green function at high energies. It is the power that 
one expects to observe in 7*7* collisions or forward-jet and Mueller-Navelet 
jet observables at ep and pp colliders respectively [22]. 

Also plotted is a second quantity labelled tOc- This corresponds to the 
position of the singularity of the gluon anomalous dimension, i.e. the power 
growth of small-x splitting functions. Though we have not really discussed 
the resummed gluon anomalous dimension, it is worth noting that ujc is 

significantly different from cOg because it contains additional corrections 

5/3 
0(as ), which arise because the effective scale for BFKL evolution in the 

anomalous dimension turns out, dynamically, to be considerably higher than 
Q^. Corrections of this form were first noticed in [31]. In general, such cor- 
rections are present for quantities involving an effective cutoff on the lowest 
accessible transverse momentum. Another example of such a quantity is the 
elastic-scattering cross section. It should be emphasised therefore that the 
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Fig. 11. Various BFKL exponents 



difference between a;^ and uJc is not an uncertainty on the BFKL exponent, 
but rather refiects differences between various processes. 

The actual uncertainty on the results can be determined by examining 
the effect of scheme changes and different approaches to the details of the 
resummation (as discussed for example in the previous section, with regards 
to the shifting of divergences) , as well as a study of solvable models [22] or 
other possible higher-order effects [29]. For ag ~ 0.2 it is about 15%. 



5. Conclusions and Outlook 

In these lectures, we have seen how to deduce many of the proper- 
ties of the BFKL pomeron. The recurrent theme has been the study of 
the collinear (and anti-collinear) limit, which gives information about the 
structure of divergences of the BFKL characteristic function at all orders. 
At NLL order the information thus obtained is sufficient to reproduce the 
true NLL corrections to a high degree of accuracy, i.e. the non-collinearly 
enhanced NLL corrections are small. Ensuring that the BFKL kernel cor- 
rectly reproduces the collinear limit at all orders leads to stable predictions 
for the high-energy power growth. The resulting resummed power is much 
more compatible with the data than either the LL or NLL values. 

For actual phenomenology, two more ingredients are required. First we 
should understand the exponentiation of the characteristic function, because 
the running of the coupling complicates the simple approach that we had at 
leading order — it turns out however that these complications are not too 
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severe [11,22,32]. 

Secondly we need to know the virtual photon impact factors, i.e. the 
coupling of a virtual photon (DIS or 7*7*) to the gluon chain. These have 
still to be worked out at NLL. When results are eventually available it is 
likely that a collinear resummation will again be needed in order to ob- 
tain stable predictions, in analogy with the situation for the characteristic 
function. 

The overall message is that despite initial fears, the large size of the NLL 
corrections BFKL is not an impediment to the use of BFKL resummation 
for predicting high-energy phenomena. But it is necessary to understand 
the origin of the large corrections, and include at all orders the physics 
which causes them. 
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